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I. INTRODUCTION 

The discovery of a discrete family of asymptotically flat particle-like solutions for 
the static, spherically symmetric Einstein- Yang-Mills (EYM) equations with SU{2) 
gauge group, made by Bartnik and McKinnon in 1988 |j|], evoked considerable interest 
in these equations and their various generalizations. The intensity of investigations 
performed in this field is shown by the latest review Q, which summarizes a decade's 
work and contains more than three hundred references to publications on the subject. 
However, there are still some problems which remain unsolved. One of the most 
interesting of these is probably the task to prove the existence of the metric oscillations 
in the origin, r = 0, neighborhood, which were found numerically during the study of 
the EYM black holes interior structure y]. The initial purpose of the present work 
was to solve this problem. Some other open questions can be found in ^. 

Let us note that there are at least two approaches to the analysis of local solutions 
of nonlinear ordinary differential equations. One of them, namely, the asymptotic 
theory of differential equations, in some cases makes a possibility to obtain a complete 
classification of solutions in a singular point neighborhood. However, the right hand 
side of the studied equation must as a rule satisfy some rather specific conditions (see, 
e.g., 1^]). Another way, known as the theory of dynamical systems, or the qualitative 
theory of differential equations, is less restrictive in this sense, though it also does not 
always lead to a comprehensive description of the solutions behavior (see, e.g., ^). 
Nevertheless, there are a number of problems in astrophysics and cosmology, which 
were solved basing on this approach (see 0, |^ and references therein) . 

In this paper, dynamical systems methods are used for the analysis of the EYM 
solutions asymptotic behavior. This enables us to prove the existence of the above 
mentioned solutions with the oscillating metric, as well as the existence of local solu- 
tions for all known formal power series expansions, and to find two new local solutions. 
Moreover, a classification of local solutions in the vicinity of the origin is obtained. 
In particular, it is shown that there exists a neighborhood of r = such that the 
metric function has a fixed sign in it. Specifically, if the limiting value of the gauge 
function equals ±1, then all real solutions belong to the Schwarzschild and Bartnik- 
McKinnon type. In other cases, the solution behavior depends on the metric function 
sign. Namely, if the metric function is positive, then all solutions possess the behavior 
of the Reissner-Nordstrom type. If, on the contrary, the metric function is negative, 
then almost all solutions are such that the metric function oscillates with its ampli- 
tude growing unboundedly as r — > 0, but the gauge function is monotonous (though its 
derivative also oscillates with the unboundedly growing amplitude). Only particular 



solutions in this case exhibit asymptotic behavior of the "anti-Reissner-Nordstroni" 
type. This result also gives the negative answer to the question stated in Q, whether 
r = is a limit point for zeros of the metric function. 

We have also considered the asymptotic behavior of solutions in the far field, r ^ 1, 
and in the vicinity of the points where the metric function tends to zero. This analysis 
leads to a discovery of two new local singular solutions and allows us to obtain some 
conclusions concerning the extendibility of solutions and their limiting behavior as the 
number of the gauge function nodes tends to infinity. 

A detailed discussion of the physical interpretation of the EYM equations solutions 
can be found in Q. 

II. THE EQUATIONS 

Recall that the space-time metric for the static, spherically symmetric EYM equations 
can be written as 

ds'^ = a'^N dt^ - N-'^ dr^ - r^ [d-d^ + sin^ -d dip^) , 

where N and a depend on r, and the Yang-Mills gauge field reads as 

A = (T2 d-d — Ti sin -d dip) w + T3 cos 1} dip, 

where Ti = -^Ti are the SU{2) group generators and Tj are the Pauli matrices, i = 1,2,3 
(see, e.g., §). 

The EYM equations in this framework take the form of two ordinary differential 
equations for the metric function N and the gauge function w. 



r^N' + (1 + 2w''^)r'^N + {l - w'^f - r"^ = 0, 

,.2\2 2,2- 



[l-uf) -r^ +r^N w' + {l-'uf)rw = {), 



(1) 



r'^Nw" 
and a decoupled equation for a: 

a' _ 2w'^ 
a r 

Since (|^) do not involve a, one can use these to obtain A'^ and w, and then solve 
the equation for a. Thus we restrict our considerations to Eqs. (|^). We also remark 
that (|^) are invariant under the transformation r — > — r; thus, in what follows we 
discuss only the region r > 0. 



For the purposes of studying the EYM solutions at finite r, it is convincing to 
rewrite (|l]) in terms of w and u = r'^N. They become 



ru 



ruw 



[l-2w' )u+ (1- 

,.2\2 



2\2 



u+il 



w 
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W 



' + {l — w ) rw = 0. 



(2) 



Recall that the only known explicit solutions of (y) are the Schwarzschild solution 



and the Reissner-Nordstrom solution 



u 



ar + r , 



w = Q, u = 1 + hr + r 



(3) 
(4) 



where a and h are arbitrary constants. 

In order to apply the theory of dynamical systems to the analysis of the EYM 
equations (0), it is necessary to write them as an autonomous system of first-order 
differential equations. Thus we introduce the function v = w' and an independent 
variable t defined by dr = rudt. After making these changes, we obtain the dynamical 
system 



ru, 



1 - 2v^) M - (1 

u+ (l -w^y 



w^?+r^ 



u, 



1 — w ) rw, 



(5) 



w = ruv. 



Notice that this system has solutions for r = and u = 0, which do not take place 
for (I). 

The first step to start analyzing (P) is to determine the critical points. It is easy 
to verify that the dynamical system (^) has the following critical sets: 



ARN^: 

RN: 

SBM^: 

W: 

RH: 



(0, -(1-^2)2, ±l,t^), 

(0,(l-u;2)2,0,«;), 

(0,0,«,±1), 

(0,0,0, w), 

[r, 0, (1 - vo^yw/l^i - w'^)'^ - r'^],w), 

(±l,0,t;,0). 

All the critical sets belong to the hyperplanes r = and/or u = 0, va which the 
conditions of the existence-uniqueness theorem do not hold for (^) . 

In what follows we shall not give a global phase portrait for (^, but we shall mainly 
concentrate on the results that have direct consequence for the EYM equations O). 



III. PRELIMINARY INVESTIGATION OF THE 
ORIGIN NEIGHBORHOOD 

Let us consider the projection of (|5[) into the hyperplane r = 0. This immediately 
leads to w = wq = const. Thus the dynamical system (|5|) reduces to 

ii = -o?u^ (l -2v^)n^, 

2 ^^) 

f = a f + liu, 

where a = \ — w^. Notice that (g) is invariant under the transformation v — > —v. 
Hence, the phase portrait will be symmetric with respect to the t(-axis. 

Since @ contains a free parameter a, it is convincing to split the analysis into 
two steps. Let us begin with a = 0. In this case, @ reads as 

u={\- 2v^\ u^ 

^ ^ (7) 

V = uv. 

One can easily solve this system. First, the critical points, which are the projection 
of the sets SBM and W , give u = Q,v = const. Next, 

C2 - Cit = - e^' - 2 / e*' dv = - e^' + i^ erffii;), and u = Civ e"^', 
V Jo V 

where Ci and C2 are arbitrary constants, and v ^ 0. Finally, for v = one has 
u = {C — t)^^, where C is an arbitrary constant, t ^ C. We remark that if n < 0, 
then the nontrivial solutions tend to zero as t ^ +cxd. In the opposite case, they tend 
to zero as t — > —00. In particular, it follows that if limj.^o^('^) = if for a solution 
of (D, then liuir^ow'ir) = 0. 

The phase portrait of (|^) is shown in Fig. Q. 

Now let us study @ for q 7^ (i.e., for wq / ±1). In this case, the system ^ has 
the following critical points: Z: (0,0), A'^: (— a^,±l), and R: (a^,0). 

The point Z is the projection of W. The eigenvalues of Z are Xu = —c? and 
\y = o? . Thus, Z is a saddle. It has four separatrices, which can be easily obtained 
explicitly. The repelling separatrices, denote them by S^ in accordance with the sign 
of V, are tangent to the eigenvector (^^ = (0, 1). They can be written as 

u = 0, v = C^ e"'*. (8) 

Here and forth the letter C, with an alphabetical subscript {Cu, C^,, etc.), denotes a 
nonzero constant. 




Figure 1: Phase portrait of the dynamical system (0). 



The attracting separatrices are tangent to the eigenvector Qu = (1)0). They take 
the form 

U = 57, v = Q. (9) 

One of these separatrices, denote it by ^i, belongs to the half-plane u < 0. For 5i, 
Cu < and t > — a~^ln|C„|. Another separatrix, 52, lies in the half-plane u > 0. It 
has Cu > and joins Z to R (see Fig. ^. 

The point R is the projection of the critical curve RN . The eigenvalues of R are 
Xu = o? and A?, = lo? . Thus, it is an unstable node. Almost all trajectories that 
approach R as t ^ — oo, are tangent to the eigenvector (u = (l, 0). The corresponding 
separatrices have the form (pi). One of them, namely, 52 joins R to Z. Another one, 
53, is defined for C^ < and t < —a~'^ln\Cu\- 

There are also two separatrices, denote them by T^, which are tangent to the 
eigenvector (^^ = (0, 1). Let us show that they have the form 



a^d 



')+oi 



(10) 



as V 
as 



0. To see this, define u., 



(n 



a }v 



'^. Now the dynamical system (P) reads 



Uy 



-a In. 



^, + 2a^v) - 2 (2a^ + Uy 



V UvV 



2c?v -|- UyV^ . 




Figure 2: Phase portrait of the dynamical system (|6|) in the vicinity of the finite 
critical points for a 7^ 0. Dashed lines represent the curves u'{v) = 0. 



For V = 0, this system has a saddle (0,0) with the eigenvalues Au„ = —a^, Xv = 2a^ 



and the eigenvectors Cu„ = (1;0), (y 



-W 



1) . The separatrices that are tangent 



to the eigenvector C,u^, belong to the line f = 0. Hence, there are no corresponding tra- 
jectories of (g). Conversely, the eigenvector C,^ determines the outgoing separatrices, 
which take the form u^ = —^a^v + o{v) as v ^ 0. This yields ([l0|) . 

Note that the same technique can be used to find the higher order terms in (10). 
This is also valid for the asymptotic solutions presented below. 

Finally, the points A represent projections of the critical curves ARN . The 
eigenvalues of A^ are A^^^ = |a^(l ± i\/l5). Thus, these critical points are repelling 
foci. It is important to note here that for all the trajectories that spiral away from A , 
the metric function u is strictly negative because of the separatrices S^, and v pre- 
serves its sign due to the separatrix Si (or, the same, because of the above mentioned 
invariance of (0) under the transformation v -^ —v). Since there are no other finite 
critical points in the half-plane n < 0, the trajectories that spiral away from the 
points A , do not have limit cycles. We remark that this can also be easily proved if 
we define the Dulac function for (g) by {u'^v)~^. Hence, u and v exhibit oscillations 



with the amphtude growing infinitely as t ^ +00. 

Fig. ^ shows the phase portrait of (^) near the points A , Z, and R. Notice 
that this portrait is drastically different from that one shown in Fig. ||. Thus, a is a 
bifurcation parameter for the dynamical system (g). 

In order to obtain the global phase portrait of (|6|), one has to study the behavior 
of its trajectories at infinity. Using the standard transform to the projective coordi- 
nates, one can find out that the system (|6|) has four critical points at the {u, v) phase 
plane boundary, namely, U^: {u = ±00, i" = 0) and V^: {u = 0,v = ±00). The 
points U^ are saddles, and V^ are saddle-nodes. The separatrices 5*3 and Si are the 
only trajectories, which approach the points U from the finite region of the phase 
plane. The points V^, besides the separatrices S^, have ingoing trajectories, which 
emanate from R. These trajectories have the form u = ^a^ f ~^ + o{v^'^) as v ^>- 00. 
The boundary of the phase plane contains two separatrices that join U^ to V and 
two separatrices that go from V^ to U~. 

IV. THE ORIGIN NEIGHBORHOOD 

A. The critical curves ARN^ 

Let us turn to the analysis of the dynamical system (^) near the critical curves ARN : 
{0,—{l — w'^y,±l,w) for u) / ±1. The excluded points also belong to the lines SBM 
and will be studied below. Notice that the curves ARN (with the points u; = ±1 
excluded) lie in the region u < 0. Hence, in a neighborhood of these curves, t growing 
to infinity corresponds to decreasing r in the EYM equations (||). 

The eigenvalues of ARN^ are A,. = -(1 - w^)^, Xu,v = ^(1 - tt'^)^(l ± za/T5), and 
A^ = 0. It should be recalled at this point that an n-dimensional critical set necessarily 
has n zero eigenvalues (see, e.g., [^). Thus, the zero eigenvalue Xw corresponds to the 
fact that ARN are one-dimensional sets of critical points. Since other eigenvalues 
have nonzero real part, ARN are hyperbolic sets. 

The eigenvalues Xu,v determine three-dimensional unstable manifolds M of the 
curves ARN , respectively. Since Xu,v are complex, the trajectories that lie on M^, 
describe oscillatory behavior of u and v. The found above trajectories that spiral 
away from the points A , are the projection of the trajectories that lie on M , in the 
plane {r = 0,w = wq j^ ±1). The separatrices 5^ and Si do also have the obvious 
counterparts for (|5|): 

r = 0, M = 0, v = Cve°\ w = wo, (11) 



and 



a2 



r = 0, u = -7, V = 0, w = wn, (12) 

respectively, where wq = const ^ ±1, and C„ < 0. These two-dimensional separatrices 
preserve the signs of u and v for the trajectories on M^. 

Recall that the trajectories that spiral away from A^ do not have limit cycles. 
Evidently, the same is valid for the trajectories on M . 

Next, due to the negative eigenvalue A^, each of the curves ARN^ has two- 
dimensional stable separatrices, which are tangent to the eigenvectors 



Ct= {l,±A{l-w^)w,-^^,±l 



w 



where the upper sign applies for ARN^ and the lower one for ARN~ . It is easy to see 
that these separatrices correspond to a one-parameter family of the EYM solutions 
that exist in the origin neighborhood. Thus, we conclude with 

Proposition 1. Let U^ be a set of solutions for the EYM equations (^) such that 
they are defined in some neighborhood of r = 0, u[r) < in this neighborhood, and 
liuij.^o'w^r) = wq < cxD, wq 7^ ±1. Then U^ is nonempty. Moreover, almost all 
solutions of (^) that belong to U^ , are monotonous for the gauge function w and 
oscillating for the metric function u. These solutions have the following properties: 

1. The amplitude of the metric function oscillations grows unboundedly as r ^ 0. 

2. The values of the metric function at the points of maximum form a sequence, 
which monotonically converges to zero as r ^ 0. 

3. The derivative of the gauge function also oscillates with the amplitude growing 
unboundedly as r ^ 0, and gets closer to zero on each cycle of the oscillations, 
but its sign remains unchanged. 

Besides these solutions, Li~ also contains a one-parameter family of local solutions of 
the "anti-Reissner-Nordstrom" type: 

u = -(1 - wlf ± 4(1 - vul)uJor + o{r), 

Wq 2 , ^ 2^ (13) 

as r ^ 0, where wq / ±1. 

Formal expansions of the form ( p!3| ) and some corresponding numerical solutions 
were found in ^. This paper was also the first to present the oscillating solutions. 
Some of their properties were analyzed in M. 



w = Wq ± r + -^^^ 27 "^ + o\r 



B. The critical curve RN 

Now let us study (|5|) in the vicinity of the curve RN: (0,(1 — w'^y,0,w) for w / 
±1. Similar to the above, the excluded points also belong to the critical sets SBM 
(and W) and will be studied below. Notice that RN (with the points it; = ±1 excluded) 
lies in the region u > 0. Hence, in a neighborhood of RN , t growing to infinity 
corresponds to increasing r in the EYM equations (^). 

The eigenvalues of RN are A^ = A^ = (1 — w^)^, A^, = 2(1 — w'^)'^, and A^ = 0. 
Therefore, all trajectories of (0) in the vicinity of RN belong to an unstable four- 
dimensional manifold; they correspond to a three-parameter family of EYM solutions. 

The eigenvalue A^, determines two-dimensional separatrices, which are tangent to 
the eigenvector (^ = (0, 0, 1, 0) and take the form 

r = 0, w = const 7^ ±1, u = a (l — |w ) + o{v ) as u ^^ 0. 

The separatrices T found above, represent their projection in the plane {r = 0,w = 
wq 7^ il)- The EYM equations (P) do not have any corresponding solution. Con- 
versely, the eigenvectors Cr = (1)0, tt;/(l — w'^),0) and Cu = (0,1,0,0) determine 
trajectories, which correspond to the EYM solutions 

u= {I - Wg)^ + uir + o{r), 

^0 2 , / 2n (14) 

as r ^ 0, where wq 7^ ±1 and ui are arbitrary constants, and the higher order terms 
contain one more parameter. 

Let us show how one can choose the third parameter in (0). The procedure will 
be similar to that one used for obtaining (|To|). Namely, consider (^) in the local 
coordinates 

u — (1 — w^)^ V 

r, Ur = , Vr = -, W. 

Then the corresponding dynamical system, which we omit for brevity, has the critical 
surface (0, Mr.,u;/(1 — uP'),w). (The other critical sets either have w = ±1 or do not 
belong to the hyperplane r = 0.) The eigenvalues of this surface are A^ = A„^ = (1 — 
nP')'^ and A„^ = A^^i = 0. It follows that all trajectories in a neighborhood of this surface 
belong to an unstable four-dimensional manifold. The nonzero eigenvalues have the 
eigenvectors C,r = (1,1 + 2i(;^,0,0) and C,vr = (0,0,1,0). Thus, all the trajectories 
assume the form 

Ur = ui + {l + 2w'^)r + o{r), Vr = —^+vir + o{r), w = wq + o{r) 



Wq 



10 



as r ^ 0, where wq ^ ±1, ui, and vi are arbitrary constants. Changing back to the 
initial variables and taking into account the above discussion, one has 

Proposition 2. Let lA^ be a set of solutions for the EYM equations (^) such that 
they are defined in some neighborhood of r = 0, u[r) > in this neighborhood, and 
liuir^o w{r) = wq < oo, Wo 7^ ±1. Then U^ is nonempty. Moreover, all solu- 
tions of (Q) that belong to U^ , form a three-parameter family of local solutions of the 
Reissner-Nordstrom type: 

M = (1 - wlf + uir + r^ + oir"^), 

w = wo + — 7y- r + w-^r^ + o[r^) 

2(1 -Wo) 

as r ^ 0, where wq, ui and W3 are arbitrary constants, wq ^ ±1. 



A formal power series expansion ( [15|) was presented in ||l^ . Some black hole solu- 
tions with this asymptotic were first found numerically in 0] . The local existence proof 
for these solutions was given in pi]. We remark that here we follow the terminology, 
introduced in Q, which is slightly different from that one used in |11| and [Q]. 



C. The critical lines SBM^ 

The lines SBM : (0,0, f,±l) are degenerate critical sets, since the eigenvalues A,., 
Am, and A^ are equal to zero. In order to study the behavior of trajectories of (^) in 
a neighborhood of these lines, we use the standard technique 0. 

First, define why w = tDibl, where the upper sign applies for SBM~^ and the lower 
one for SBM~ . Now the lines SBM are transformed to the v-axis. Next, introduce 
the local coordinates 

ru = -, u, V, Wu = — , (16) 

u u 

in which the dynamical system (|5|) reads as 

fu = {2v^+K)ru, 

u = (l-2v^ -K)u, 

(17) 

V = {1 + K)v + {1± uwu){2 ± uWu)ruUWu, 

Wu = -{I- 2v'^ - K)wu + ruV, 

where K = [{2 i: uwu)'^w'^ — r'^]u, and an overdot stands for derivatives with respect 
to r defined by dr = udt (thus, r = Inr + const). 



11 



The system (O) has one critical set in the hyperplane n = 0, namely, the r„-axis, 
which is an unstable hyperbolic line. The corresponding eigenvalues are A^^ = 0, 
Am = Ai, = 1, and A„,^ = — 1. The two-dimensional ingoing separatrices 

r„ = const, u = 0, v = 0, Wu = C^ e""^, 

which are tangent to the eigenvector C«)„ = (0, 0, 0, 1), belong to the hyperplane u = 0. 
Hence, they do not correspond to any trajectories of @. In their turn, the eigenval- 
ues Xu and A^, which have the eigenvectors Cm = ("'^mj 1)0,0) and (^ = (0)0,2,rM), 
determine the outgoing three-dimensional separatrices 

ru = ro -rQU + o{u), v = 2viu + o{u), Wu = vqVi u + o{u) 

as n ^ 0, where tq and vi are arbitrary constants. This implies 

Proposition 3. All solutions of the EYM equations (0) such that \mir^ow{r) = ±1, 
belong to a two-parameter family of local solutions of the Schwarzschild type: 

u = uir + r + o(r ), 

(18) 
w = ±1 + W2r + o{r ) 

as r ^ 0, where ui and W2 are arbitrary constants. 

In particular, these local solutions describe the behavior of the Bartnik-McKinnon 
particle-like solutions [||] (for ui = 0) and the black hole solutions of the Schwarzschild 
type [^ in the vicinity of the origin. 

It is interesting to note that the separatrices that are tangent to the eigenvector Cm, 
can be written as 

1 + uiru 



' II 



0, Wu = 0, 



where ui is an arbitrary constant, and Tm 7^ 0. Obviously, these separatrices corre- 
spond to the Schwarzschild solution (^. 

We also remark that the coordinates (|l^) enable us to resolve the degeneracy of 
the T;-axis along the u-direction. Analysis of the r- and tu-directions leads to the same 
conclusion for the EYM equations as stated in Proposition 0. 

It is necessary to underline here that we discuss only real EYM solutions, though 
the EYM equations also possess complex solutions. For example, a study of (^) in 
the local coordinates {r,u/r'^,v,w/r), leads to a discovery of complex EYM solutions 
of the form 

u= - {1 + iwj) (r^ + WQWir^^) + o{r^), 
ty = iuq + Wir — ^wor'^ + o{r^) 

12 



as r ^ 0, where wq = ±1 and wi = ib(-v/3 it i-v/5)/4- These solutions do not have free 
parameters. 

D. The critical line W 

The eigenvalues of the critical line W: {0,0,0, w) are A^ = 0, A^ = —(1 — w'^)'^, 
A^ = (1 — w^)^, and A^ = 0. Hence, W is degenerate for any w. The eigenvalues Xu 
and Xy are nonzero and have different signs whenever w ^ ±1. In this case, W is 
an unstable critical set with the outgoing two-dimensional separatrices ([ll|) and the 
ingoing two-dimensional separatrices (12). Investigation of (||) in the vicinity of W 



gives the same result for the EYM equations as already stated in Proposition |3|. By 
this reason we omit the discussion. 

Thus, we have obtained a description of the EYM solutions that have finite values 
of the gauge function w in the origin neighborhood. It follows from our considerations 
that for all these solutions limr^o w{r) = wq < oo. A standard analysis of the 
dynamical system (|5|) in the projective coordinates {rz, uz, vz, z = l/w) reveals that 
Eqs. (y) do not have solutions such that limr^o^('") = oo. Hence, the results of 
this section can be summarized in the following classification of the EYM solutions, 
defined in the vicinity of the origin. 

Theorem. All real solutions of the EYM equations ^, defined in a neighborhood of 
r = 0, belong to one of the following disjoint classes: 

1. Wq = ^1. In this case, all solutions are of the Schwarzschild and Bartnik- 
McKinnon type (ll^ 



2. Wq 7^ ±1, and the metric function u is negative in some neighborhood of r = 0. 
In this case, almost all solutions are such that the metric function oscillates 
with the unboundedly growing amplitude as r — > 0, but the gauge function is 
monotonous (though its derivative also oscillates with the amplitude growing 
infinitely). Only particular solutions in this case exhibit asymptotic behavior of 
the "anti-Reissner-Nordstrom" type (^ 



3. Wq 7^ ±1, and the metric function u is positive in some neighborhood of r = 0. 
In this case, all solutions belong to the Reissner-Nordstrom type (15). 



We remark that this classification scheme explains why almost all interior black 
hole solutions, found numerically in 0], exhibit oscillatory behavior of the metric. 
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V. CRITICAL SETS FOR r 7^ 

A. The critical surface RH 

Now let us discuss the remaining critical sets of (^). The surface RH: (r, 0, (1 — 
w'^)rw/[{\ — vlP')'^ — r^], w) is an unstable hyperbolic set. The eigenvalues of RH are 
\r = 0, A„ = -[(l-w^)^-r^], At, = (l-w^)^-r^, and A^ = 0. The three-dimensional 
separatrices that are tangent to the eigenvector Qy = (0,0, 1,0) read as 

r = ro = const, n = 0, f = — + C„ e'^*, w = wq = const, 

where /3 = (1 — Wq)^ — Tq 7^ and 7 = (1 — w^)rQWQ. Obviously, they do not 
correspond to any EYM solution. In their turn, the three-dimensional separatrices 
that are tangent to the eigenvector 

_ / G {2F'^ - [F^ + (1 - 3w;2)r2]r2}Fw Frw\ 

where F = 1 — w^ and G = F'^ — r'^, correspond to a two-parameter family of local 
EYM solutions. It is convenient to fix one of these parameters and to write down 
these solutions as follows. 

Proposition 4. For any fixed Vh > 0, the EYM equations (§) possess a one-parameter 
family of local solutions of the form 

{l-wlf-rl 
u = s + o(s), 

" = "^+ :i-^Jplrg ^ (19) 



(1 - w%)rhWh 

,,2^/0^ _ .,,2^4 _ rn _ ,„2\3 , n _ q,„2\^2l^2l 



(1 - wi){2{l - wif - [(1 - wiY + (1 - ^wpriVDwu ^2 ^ ^(^2) 



4[(1 - wlf - rl 



213 



as s = r — r/i — > 0, where Wh is a constant, satisfying \1 — wf^\ ^ r^. 



The local solutions (19) represent the EYM solutions in the vicinity of a regular 



horizon. For the black hole solutions, this is either an event horizon (if r/j > |1 — w'^ 



or an interior Cauchy horizon (if r/^ < |1 — wf^\). The first existence proof for these 
local solutions was given in [jl2|. Some black hole solutions with an interior horizon 
were found numerically in 0. 

Note that RH transforms to the line W for r = 0. 



14 



B. The critical lines DH^ 

The lines DH : (±1,0,^,0) are degenerate. All their eigenvalues are equal to zero. 

Since the EYM equations @ are invariant under the transformation r — > — r, we shall 

study (P) only in the vicinity of the line DH^ . 

It is convenient to use the local coordinates 

u w 

r = r-l, Ur = ^, V, Wr = —, 

in which the dynamical system (|5|) reads as 

f = (1 + f) fUr, 

iir = \2- (2 + r + 2rv^) Ur + (l + 2w^) r - f^wf\ Ur, 

V + Wr — UrV + 2vw^) r + r Wj. + [1 + VWrjr w^, 
Wr = {I + f){v — Wr) Ur, 

where an overdot stands for derivatives with respect to r defined by dr = f dt. 

There are two critical sets of ( po|) in the hyperplane f = 0, namely, the point 
D: (0,1,0,0) and the line L: {0,0,— ^Wr,Wr). The eigenvalues of D are Xf = 1, 
A^^ = —2, and At,^„,^ = — ^(3 ± i-v/3). Thus, D is a saddle. The outgoing one- 
dimensional separatrices 

r = --— "7^ -, Ur = l, V = 0, Wr = 0, 

1 + Ur e ^ 

which are tangent to the eigenvector Cf = (1)0,0,0), correspond to the extreme 
Reissner-Nordstrom solution 

w = 0, u = {l-rf. (21) 

Besides this, in the vicinity of D there exists a stable three-dimensional manifold. 
In order to figure out whether this manifold belongs to the hyperplane f = 0, one 
may study a projection of ( [20| ) into f = const. It occurs that the critical point 
(1,0,0) exists for any f and has the eigenvalues A„^ = — 2 — f and A^^^,^ = — ^ (3 + 
f) lb iy/{3 + f)(l + 3f). Hence, for any f > —1/3 all the eigenvalues have negative 
real part, and A^, and A^;^ are complex conjugate. Thus, the stable three-dimensional 
manifold does not belong to the hyperplane f = 0, and the trajectories on this manifold 
correspond to a two-parameter family of local EYM solutions. 



It is interesting to note that the projection of (20) into the plane {r = 0,Ur = 1) 
gives a system of two linear differential equations 

V = —2v — Wr, 

(22) 

Wr = V — Wr, 
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which can be easily solved: 



Wr 



Cx cos ('^r) - ^(Ci + 2C2) sin {^t 
C2 cos (^r) + ^{2Ci + C2) sin (^r 



exp (-|t) 
exp (-|r) 



where Ci and C2 are the constants of integration. Thus, the projection of (pQ) into 
(f = 0, Ur = 1) represents linear oscillations of v and 11;^ with infinitely many zeros. 

It is easy to see that f — > corresponds to r — > —00, so that w' diverges as r — > 1. 
Conversely, r goes away from 1 as r — > +00, and these solutions tend to the extreme 
Reissner-Nordstrom solution. Thus, we have 

Proposition 5. There exists a neighborhood of r = 1, in which the EYM equa- 
tions (0) have a two-parameter family of local solutions with the gauge function oscil- 
lating with infinitely many zeros, and the metric function tending to zero. 

These solutions can be treated as a description of the limiting behavior of the EYM 
solutions in the vicinity of r = 1 as the number of the gauge function nodes tends to 
infinity. The limiting behavior of the EYM solutions was first studied in [13, 14, 15]. 



Solutions that exhibit oscillations of w were first discussed in [^. But, in addition to 
the results of []lq ], we see that the gauge function may have infinitely many zeros not 
only to the left of r = 1, but also to the right. 

Finally, the line L has the eigenvalues Af = 0, A^^ = 2, A^, = —2, and Xw^ = 0. 
Thus, L is a degenerate set. The eigenvectors C^^ = {0,1, j^Wr,—jWr) and (y = 
(0,0, 1,0) determine two-dimensional separatrices, which lie in the hyperplane f = 0. 
Thus, they do not correspond to any EYM solution. Further investigation of the line L 
did not reveal any trajectories of (^) that have corresponding EYM solutions different 
from the discussed above. 

VI. SOLUTIONS WITH A SINGULAR HORIZON 

A typical EYM solution cannot be continued to the far field, since it has a singular 
horizon, i.e., a point, at which the metric function tends to zero, the gauge functions 



stays finite, but its derivative diverges [17, 12]. This fact was first noticed in |13], 



and a power series expansion, describing the behavior of the EYM solutions in the 
vicinity of a singular horizon was given. Let us show how one can obtain singular 
EYM solutions basing on dynamical systems methods. This will also lead us to a 
discovery of two new local solutions. 
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Let us rewrite the dynamical system (pi) as 



r = ruz 
u = 

z = 



2\2 ^2 „2^2 



(2-z')u+{l-wYz^-r^z 



u, 



+ [l — w j — r — (l — w) rzw 



(23) 



w = ruz, 



where z = 1/v, and an overdot stands for derivatives with respect to r defined by 
dt = z^ dr. 

For 2 = 0, the critical points of (P3| ) form a degenerate plane (r, 0,0, if). In 
order to study ( [2^ ) in the vicinity of this plane, we introduce the local coordinates 
{r,Uz = u/z'^, z,w), in which (B^) can be written as 



r = ru^-z 



2\2 



{2 - 3z^) u^ - {l - w^) +r^ + 2{l-w^) 



rzw 



UzZ 



2\2 



+ (1 - W^) 



(1--^) 



rztM 



(24) 



w = ru,z, 



where an overdot stands for derivatives with respect to t. 

The system ( P4|) has two critical sets in the hyperplane z = 0, namely, the 
plane RWi: (r, 0,0, w) and the surface SH: (r, [(1 — w^)'^ — r^]/2,0,ti;), which are 
nondegenerate whenever 



r ^\1 — w \. 



(25) 



The eigenvalues of RWi are A,. = 0, A^^ = (1 — iiP'Y ~ '"^j ^z — ~^uzj and A^ = 0. 
Thus, RWi is an unstable hyperbolic set. One can easily see that if (|25| ) holds, then 
the eigenvectors Cu^ = (0, 1,0,0) and Cz = (0,0, 1,0) determine the three-dimensional 
separatrices 



r = ro = const, Uz 



2 + Cue-Pi' 



z = 0, w = wq = const, 



and 



r = ro, Uz = 0, z 



P 



J + Cz e/3* 



, W = Wo, 



respectively, where /3 = (1 — Wq)"^ — Tq / and 7 = (1 — ^0)^0^0- Clearly, these 
separatrices do not correspond to any EYM solution. 
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Next, the eigenvalues of SH are A,. = 0, A^^ = Xz = —[(1 — w'^)'^ — r^], and 
Xw = 0. Hence, if ( p5| ) holds, then all trajectories of ( |2^ ) in the vicinity of SH belong 
to a four-dimensional manifold. The eigenvalues Xu^ and A^ have the eigenvectors 
Cmz = (0)1)0)0) ^-iid (^2 = (0,0,1,— ^r). Thus, all trajectories in the vicinity of SH 
take the form 

r = ro + o{z), Uz = - + uiz + o{z), w = wq- —z + o{z) 

as z ^^ 0, where (3^0, and ui is an arbitrary constant. This leads to 

Proposition 6. Let vq > be a point such that 

lim u{r) = 0, lim w{r) = wq < oo, lim w (r) = co, 

r— >ro r^ro r^ro 

and (3 = {1 — Wq)'^ — Tq / 0. Then all solutions of the EYM equations @ in the 
vicinity of rg have the form 

2(3 
u = — s^ + uis^ + o(s^), w = wq rt y/vQ s + o{s) (26) 



as s = \JrQ — r -^ 0, where ui is an arbitrary constant. These solutions do not have 
other parameters besides wq and ui . 

Thus, local solutions (|26|) exist in the vicinity of any point ro > 0, ro j^ \1 — Wq\. 
Notice that solutions (|l9|) were determined by the trajectories that belong to a 



three-dimensional manifold. Unlike them, singular solutions (26) correspond to the 
trajectories that lie on a four-dimensional manifold. It follows immediately that in the 
vicinity of an arbitrary point tq > such that lim^^j-Q u{r) = 0, lim^^ro w^r) = wq, 
and rQ ^ \1 — Wq\, almost all solutions of the EYM equations (|2|) exhibit asymptotic 
behavior (p6|) and, therefore, cannot be continued toward r = oo. 

Let us also mention that it follows from (^) and the above analysis of the critical 
sets RH and DH that almost all EYM solutions, defined at an arbitrary finite point 
ro > 0, can be continued to the left for all r < tq. A detailed investigation of 
extendibility of solutions of the EYM equations can be found in B. 

Now let us study (23) in the vicinity of the curve r = |1 — i(;^|. We start with 
\uj\ < 1. In the local coordinates 

1 



rz = -{r-l + w), Uz = ^, z, w, 
z z-^ 
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the dynamical system (p^) reads as 



(l-w^) {z + 2w)uz+ 2{z + w)u;,z - [l ~ w^) 



2\2 

w w 



1 



w 



[2 + zw)ri 



,,2\2 



'1 



2 (1 - 2z^) n^ + 3 (1 - w^) w+{l- w^) (4 + 3zw)r^ + 2 

w ) (2 + zt(;)r^ + r^z 



W ) W — UzZ + (l 



r^z 



Uz 



(27) 



It; 



(l - -u;^ + r^z) UzZ, 



where an overdot stands for derivatives with respect to r defined hy dr = z dt. 

The system (^) has six critical sets in the hyperplane z = 0, namely, LW: 
{0,0, 0, w), LR^: (r^,0,0,±l), CR: {-^{l - w^)w, 0,0, w), RUi: {-{1 - w'^)w,^{l - 
w'^)'^w,0,w), and RU2: (-|(1 -w'^)w,^{l - w'^)'^w,0,w). Analysis of the first four 
critical sets does not reveal any trajectories of (^) that have corresponding EYM 
solutions. Thus we discuss only the curves RUi and RU2- 

The eigenvalues of RUi are A^^ = Xu^ = (1 — w'^)'^w, Xz = —(1 — uP')'^w, and 
A^ = 0. Thus, RUi is an unstable hyperbolic critical set whenever 



w^O,±l. 



(28) 



In this case, all trajectories on a three-dimensional manifold, determined by Xr^ 
and Xu^, are tangent to the eigenvector Cr^,u^ = (Ij— (1 — ^^);0,0), which defines 
the two-dimensional separatrices 



Uz 



-^a Wo — avz 



z = 0, w = wq ^ 0, ±1. 



One of these separatrices joins RUi to CR. However, the hole manifold belongs to 
the hyperplane z = 0. Thus, the trajectories on it do not correspond to any EYM 
solution. 

Unlike this, the two-dimensional separatrices that are tangent to the eigenvector 

Qz= {\ [l-w') {3-^w') ,-\ [l-w^f {2-3w^) ,l,-\ [l-w')) , 

take the form 



Tz = —awQ -|- I a (3 — St^o) z + o{z), 
Uz = \ a^WQ — I a^ (2 — 3wq) z + o{z), 
w = wq — ^2 OLZ + o{z) 

as z — > 0, where wq / 0, ±1, and thus have corresponding EYM solutions. 



(29) 



19 



Finally, the eigenvalues of RU2 are Xr^ = —{I — w'^)'^w, Xu^ = 3(1 — nP')'^w, 
Xz = —2(1 — w'^)'^w, and A^ = 0. Hence, RU2 is also an unstable hyperbolic critical 
set whenever (|2q ) holds. The eigenvector (^„^ = (1,— (1 — w^),0,0) defines the two- 
dimensional separatrices 

Uz = —OLTz-, 2; = 0, u; = -Wo 7^ 0, ±1. 

One of them joins BXJ2 to LW . Besides these separatrices, there is also a three- 
dimensional manifold, defined by the eigenvalues Xr^ and A2. Almost all trajectories 
on this manifold are tangent to the eigenvector Cr^ = (1) — 3(1 — i(;^),0, 0), which 
determines the separatrices 

Uz = \oL Wo — 3as -|- o{s) as s = r^ -|- | awo — > 0, z = 0, w = wq, 

where wq 7^ 0,ibl. In addition, there exist two-dimensional separatrices, which are 
tangent to the eigenvector 

C. = (I) (1 - ^') (13 - 12u;2) , -|j (1 - n^'f (11 - 9^') , 1, -! (1 - w^ 
and take the form 

-| avjQ -I- Jj a (13 - 12^0) z + o{z), 



r 



z 



U,_ 



a^ WQ -^a^ [11 -9w'q)z + o{z), (30) 



2 " "^u 40 



w = wq — '^az + o{z) 

as 2: — > 0, where wq 7^ 0, ±1. These separatrices, together with (|29|), have correspond- 
ing singular EYM solutions. Conversely, all trajectories, defined by the eigenvalues A^^ 
and Xr^ , belong to the hyperplane z = and do not have counterparts neither for the 
dynamical system (p3|), nor for the EYM equations. 

Analysis of ( P3|) in the vicinity of the curve r = — (1 — u;^) for lii;] > 1 is completely 
analogous to the previous case. The dynamical system (^), written in the local 
coordinates 

rz = -{r + l-w'^), Uz, z, w, 

z 

has the same critical sets in the hyperplane z = 0, as (]27|), to the exclusion of the 
curves RUi and RU2, which in this case have the opposite sign of Uz- Asymptotic 
formulas (^) and (^) convert to 



Tz = —awo — I a (3 — 5wl) z + o{z), 
Uz = —^a^wo — \a^ [2 — 2,Wq) z + o(z). 



w = wq + ^az + o{z), 



20 



and 



^z - 


= -\awQ - 


|5"(13- 


\2wl) z 


+ < 


AA. 


Uz - 


^-\o?WQ- 


-|l"Mll 


- 9«^o) 


2 + 


o{z), 


W - 


= -Wo + f a^; 


+ o(z) 









as z — > 0, respectively. Recall that a = 1 



-Wn 



-ro 7^ 0,1 here. Combining these 



solutions with (P9|) and (pO|), we get 



Proposition 7. /n the vicinity of any point rg > 0, vq ^ 1, the EYM equations 
have solutions of the form 



and 



u 



u 



1 c\ 

±—iwoWi2S^ + 0{s^), 



w 



Wo ±^/r^s + o{s), 



Wo ±Wi2S + 0{s), 



as s = \fro^^ —>■ 0, where ro = \1 — Wq\, ^ 
These solutions do not have free parameters. 



sgn(l - wl), and wu = \/3ro/2. 



To the best of the author's knowledge, the presented local solutions are new. 

Notice that Proposition excludes the cases tq = 1 and ro = 0, which correspond to 
wo = 0, ±1. Analysis of the first one reveals complex solutions of the EYM equations. 
We do not discuss them here, since their physical interpretation is unclear. The latter 



case has already been discussed in Sec. IV, 



VII. SOLUTIONS IN THE FAR FIELD 

The behavior of the EYM solutions in the far field, r ^ 1, was studied in great details, 
see 1 18] and references therein. In this section, we briefly give another existence proof 
for the asymptotically flat solutions and obtain a description of the limiting behavior 
of the EYM solutions as the number of the gauge function nodes tends to infinity. To 
implement this task, we return to the EYM equations (|l]), but we change r to z = 1/r. 
Next, we rewrite (|^) as a dynamical system of the form 



(31) 



i = z'^N, 






N = 


'-l+^l + 2z\^)N+{l-w^)\^' 


zN, 


V = 


1-3N- {l-w'^^z'^ 


zv — [l — w )w 


W = X 


z'^vN, 
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where v = w'{z), and an overdot stands for derivatives with respect to t defined by 
dz = z^Ndt. 

The dynamical system ( |3l|) has two critical sets in the hyperplane z = 0, namely, 
the planes AF : (0,A^, ^,±1) and OS: (0,A^, f,0). Both of them are degenerate. 
Thus we perform the standard procedure of their investigation for finite N and v. 

A. The critical planes AF"^ 

In this case, we introduce the local coordinates (z, N, v, Wz = w/z), where w = w ±1; 
here the upper sign applies for AF~^ and the lower one for AF" . Now ( |3l|) can be 
written as 

z = zN, 

AT = [-1 + (1 + 2z^f2) N + {2± zw^fzW,] N, 

v= [1-3N -{2± zwzfz'^wl] v + {l± zw^){2 ± zwz)wz, 

Wz = {v- Wz)N, 

where an overdot stands for derivatives with respect to r defined hy dr = z dt. 

The dynamical system ( p^ has two critical lines in the hyperplane z = 0, namely, 
Zi: {0,1, WzjWz) and Z2: {0,0,—2wz,Wz)- The eigenvalues of Zi are Xz = \n = Ij 
A^ = —3, and A^^ = 0. The eigenvector Cu = (0,0,-2,1) defines the ingoing two- 
dimensional separatrices 

z = 0, N = I, Wz = Wq — |f, 

where wq is an arbitrary constant. Since these separatrices belong to the hyperplane 
z = 0, they do not correspond to any trajectories of (p|). 

The eigenvalues A^ and Aat have the eigenvectors Cz = (l^O, ±|w^, ib|u)^) and 
Cat = (0,1, —^Wz, —^Wz) , where the signs in C,z correspond to the signs in the right 
hand sides of (^2[). Thus, the three-dimensional outgoing separatrices take the form 

A^ = l + niz + o{z), 
V = wi + |(±t(;i — ni)wiz + o{z), 
Wz = wi + j{±wi - ni)wiz + o{z) 

as z — > 0, where rii and wi are arbitrary constants. Clearly, these separatrices corre- 
spond to a two-parameter family of the asymptotically flat solutions of (|l|). 

Proposition 8. The EYM equations (|l|) possess a two-parameter family of solutions 
such that limr_»oo w{r) = Woo = ±1. All these solutions have the form 

N = l + n^ir~^ + o{r~^), 

4^ 



w = Woo + W-ir ^ + j{wooW^i - n^i)w-ir "^ + o{r ^) 
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OS r — > oo, where n_i and W-i are arbitrary constants. 

Next, the eigenvalues of Z2 are Xz = 0, \n = —1, A^, = 1, and A^,^ = 0. Hence, Z2 is 
an unstable degenerate set. The eigenvalues (n = (0, 1, — Gtt;^, Stt;^) and Cn = (0, 0, 1,0) 
determine two-dimensional separatrices, which belong to the hyperplane z = 0. Thus, 
they have no corresponding trajectories of (|3l|). Closer analysis of Z2 did not reveal 
any trajectories of (|3^) that correspond to EYM solutions. 



B. The critical plane OS 



In this case, we study (|3l| ) in the local coordinates {z,N,v,Wz = w/z), in which (|3lD 
may be written as 



z = zN, 

N = [-1 + {I + 2z'^v^) N + {1 - z^w^f z^]N, 
{,= [l-3iV- {I - z^w^zf z^]v - {l-z^wl)wz 

Wz = {v - Wz)N, 



(33) 



where an overdot stands for derivatives with respect to r. 

The system ( |33| ) has two critical sets in the hyperplane 2 = 0, namely, the point 
P: (0,1,0,0) and the line Z3: (0, 0, ii;^, w^). The eigenvalues of P are Xz = Xn = 1 
and Xv^w:, = ~2('^ ^ «\/3). Thus, P is a saddle. The eigenvectors Qz = (1,0,0,0) and 
Cat = (0, 1,0,0) determine the outgoing two-dimensional separatrices 

N = l + niz + z'^, w = 0, Wz = Q, 

where ni is an arbitrary constant. Obviously, these separatrices correspond to the 
Reissner-Nordstrom solution (^). 

Next, trajectories that belong to a stable two-dimensional manifold defined by the 
eigenvalues A^ and A^^, spiral toward P as r ^ -|-oo. These solutions may be written 
down explicitly, since for z = and A^ = 1 the system ( |33|) reads exactly as (|22|) 
with Wr replaced by Wz- However, the whole manifold belongs to the hyperplane 
z = 0, so that the trajectories on it do not correspond to any EYM solution. One may 
treat these trajectories as a description of the limiting behavior of the EYM solutions 
as the number of the gauge function nodes tends to infinity. 

Finally, the line Z3 is an unstable degenerate critical set. The eigenvalues of Z3 
are A^ = 0, Aat = —1, A^, = 1, and A^^ = 0. One can easily see that the eigenvectors 
C,N = (0, 1, |tL'2,0) and C,y = (0,0,1,0) define two-dimensional separatrices, which 
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belong to the hyperplane z = 0. Thus, they do not have corresponding trajectories 



of (31). Additional study of Z3 did not reveal any trajectories of p^ ) that have 



corresponding EYM solutions. 

Let us mention in conclusion that though our investigation was restricted to local 
solutions of the EYM equations, dynamical systems methods can also be used for the 
analysis of the solutions global behavior. This will be the subject of a forthcoming 
publication. 
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